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Abstract 

Nonperturbative QCD contributions to the inclusive semileptonic 
decay of the B meson consist of the dynamic and kinematic compo- 
nents. We calculate the decay width in an approach based on the 
light-cone expansion and the heavy quark effective theory, which is 
able to include both components of nonperturbative QCD contribu- 
tions. The kinematic component results in the phase-space extension 
and is shown to be quantitatively crucial, which could increase the 
decay width significantly. We find that the semileptonic decay width 
is enhanced by long-distance strong interactions by -|-(9 it 6)%. This 
analysis is used to determine the CKM matrix element \Vcb\ with a 
controlled theoretical error. Implications of the phase-space effects 
for the nonleptonic decay widths of 5-hadrons are briefly discussed. 
The experimental evidence for the phase-space effects is pointed out. 
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1 Introduction 



A direct goal of studying the inclusive semileptonic B meson decay B — >■ 
ei'e^ is to determine the standard model parameter |V"c6| accurately. The 
semileptonic decay width can be expressed as 



lc\Vcbf + 'yu\Vubf ■ (1) 



SL 



The first term in results from the b ^ c transition. The second term in 
(|l]) is due to the h ^ u transition and is negligible in comparison with the 
first term since 7^ ~ 7c and \Vub\ <^ \Vcb\- The semileptonic decay width is 
determined by two measured quantities: the inclusive semileptonic B decay 
branching ratio Bsl and the B meson lifetime tb, 

r.. = ^. (2) 

Therefore, the CKM matrix element \Vcb\ can be determined through 

IT. |2 ^SL Bsl /„n 

\^cb\ — — , [O) 

Ic IcTb 

with the theoretical input 7c. 

Theory is needed to calculate 7c and to understand quantitatively uncer- 
tainties in this calculation. The main obstacle to this end is the difficulty 
of taking into account nonperturbative QCD effects on the underlying weak 
decay process. 

In recent years a heavy quark expansion approach to inclusive B decays 
has been developed [1-7] to account for nonperturbative QCD effects. This 
approach is based on the operator product expansion and the heavy quark 
effective theory (HQET). An operator product expansion on the time ordered 
product of two currents is performed. The momentum of the incoming b 
quark is written as pb = ^bV + k {mi, stands for the b quark mass and v 
the B hadron velocity) and the residual momentum, k, is expanded in. For 
keeping track of the rub dependence of matrix elements, the b quark operators 
in full QCD are matched onto those in the HQET. The leading term of the 
expansion coincides with the free quark decay model. The next terms are 
computed in powers of I /nib, where no l/rrih term appears. 

The calculations P, ^ in the heavy quark expansion approach claimed 
that nonperturbative QCD contributions decrease the semileptonic decay 
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width by a few percent with respect to the free quark decay width. There are, 
however, theoretical hmitations in this approach [Q. The operator product 
expansion breaks down for low-mass final hadronic states. In particular, the 
endpoint singularities of the lepton spectra indicate a failure of the operator 
product expansion. Moreover, the truncation of the expansion enforces the 
use of quark kinematics rather than physical hadron kinematics. Describing 
the lepton spectra demands a resummation of the heavy quark expansion 
1^. There remains a need to clarify the consequence of the theoretical lim- 



itations for the calculation of the semileptonic decay width, as it is desirable 
to improve the accuracy in the independent determination^ of \Vcb\ from the 
inclusive semileptonic B decay with theoretical refinements. 

The resummation of the heavy quark expansion introduces |TD[ a distri- 



bution function ("shape function") of the b quark in the B meson, which 
incorporates nonperturbative QCD effects. A similar distribution function 
arises |Tl|, |I2[ from the light-cone dominance in the inclusive semileptonic B 
meson decays. The introduction of the distribution function eliminates the 
theoretical difficulties mentioned above, namely the endpoint singularities 
are absent and the use of physical hadron kinematics is allowed (but does 
not arise "for free" ) . 

In this paper we will use the light-cone approach [Tl|, |12| to calculate the 



nonperturbative QCD contributions. This approach describes the decay by 
using the light-cone expansion and the HQET, which provides a theoretical 
justification for the DIS-like parton model The predicted electron energy 



spectrum agrees well with the experimental measurement [1^. The use of 
physical hadron kinematics is built into this approach explicitly, so that both 
dynamic and kinematic effects of nonperturbative QCD are properly taken 
into account. The latter is shown to be quantitatively crucial. We find an 
about 9% enhancement of the semileptonic decay width with respect to the 
free quark decay width by nonperturbative QCD contributions, in contrast 
to the results obtained in the heavy quark expansion approach. 

The reason of the enhancement is the following. There are two com- 
ponents - dynamics and kinematics - of nonperturbative QCD effects on 
inclusive semileptonic B decays. First, the decay dynamics deviates from the 
free quark decay dynamics as quarks are confined in hadrons and can never be 



TUsing different experimental and theoretical methods, \Vcb\ can be also determined 
independently from exclusive semileptonic B decays. 
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free. However, the dynamic deviation changes the decay width only shghtly 
since the b quark inside the B meson is almost on shell. Second, the decay 
kinematics gets changed. The phase space extends from the quark level to 
the hadron level (the detailed formulas will be given below in (p!^)-([T5|) for 
hadron kinematics and in (|16D and (|l^) for quark kinematics), shown in Fig.|l] 
for the 6 — c decay. The phase-space extension arises from the difference 
in the B meson and b quark masses and the fact that the mass of the decay 
product quark is fixed in the free quark decay picture, while the mass of the 
final hadronic state is actually changeable. The phase-space effect is a dom- 
inating factor, as indicated by the replacement of the 6-quark mass with the 
S-meson mass in the decay rate. It is thus important to include this type of 
contributions to the decay width. Consequently, it is conceivable that the net 
effect of nonperturbative QCD enhances the semileptonic decay width. The 
negative contribution found in previous studies in the heavy quark expan- 
sion approach is just a reflection of incompleteness of the calculation, which 
fails to take into account, in particular, a large part of nonperturbative QCD 
contributions due to the phase-space effect. 

We will describe the approach in section 2 and analyse nonperturbative 
QCD contributions in comparison with the heavy quark expansion and ex- 
tract then \Vcb\ from the inclusive semileptonic B meson decay in section 3 
and finally conclude and discuss in section 4. 



2 Approach 

The semileptonic decay width can be split into two parts: one, denoted by 
^nonpert, includcs uouperturbative QCD contributions, the other results from 
perturbative QCD corrections to the decay width, denoted by Tp^rt- Namely, 

^nonpert ~l~ Tpert- (4) 

Nonperturbative QCD effects are contained in the hadronic tensor Wf^^- 
It can be written in terms of a current commutator taken between B states: 

W,. = f d'ye''^-y{B [j,{y),jlm B), (5) 

where q is the momentum transfer to the final lepton pair. | B) refers to 
the B-meson state with energy Eb and is normalized according to {B\B) = 
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It is well known that integrals like the one in eq.(5) are dominated by 
distances where 

< < 4- (6) 

For inclusive semileptonic B-meson decays, the momentum transfer squared 
is timelike and varies in the physical range 

0<g^<(M5-Mx_)^ (7) 

where Mb and Mx^i„ represent the B-meson mass and the minimum value 
of the invariant mass of the hadronic final state, respectively. Due to the 
large B-meson mass, extended regions of phase space involve large values of 
(see also Fig.l). Therefore, the decay is dominated by light-cone distances 
between the two currents in eq.(5). This allows to replace the commutator 
of the two currents with its singularity on the light cone times an operator 
bilocal in the b quark fields. Furthermore, the light-cone dominance enables 
us to expand the matrix element of the bilocal operator between B-meson 
states in powers of Aq^b/q"^- The leading nonperturbative effect is described 



by a distribution function |IT|, |12 



/(O = / d{y.PB)e^^'''-{B\b{OWB{l-l,)b{y)\B) |,.=o, (8) 

where Pb denotes the four-momentum of the B meson. /(^) is the probability 
of finding a 6-quark with momentum C,Pb inside the B meson. The hadronic 
tensor can be expressed in terms of the distribution function: 

W,, = A{S,^,p-ie,^,p) I dU{Oe{iPB,-qo)mPB-q?-K]{iPB-qTP^. 

(9) 

where rric is the charm quark mass. 

^nonpert is Calculated in this approach by integrating the differential decay 
rate in the B rest frame, 



r_,e. = jdE,j di I "1^^ (10) 



d^T 

dq^ I dMt- — - „ . 

^ J ^dE^dq^dMl' 

with the differential decay rate 
dE^dqHMl Stt^b f^^y^^^\ Mb) ^ 
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where Ee is the electron energy and we have neglected the electron mass. 
Mx denotes the invariant mass of the final hadronic state. The dimensionless 
variables ^± reads 



= ^°^Vq^ + "^-. (12) 
Mb ^ ^ 

Note that there appears in the differential decay rate (|lTD the B meson mass 
rather than the b quark mass. The integration limits are specified by hadron 
kinematics: 

Mb( Ml \ 
0< Ee <^(l-^^), (13) 

(14) 

Ml_< Ml <{MB-2Ee)(Ms-^^, (15) 

which define the hadron level phase space shown in Fig.l. For comparison, 
we also write down here the kinematic boundaries for the free quark decay: 

0< Be <^fl-4l (16) 
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0< <2EJm, ^ , (17) 

which are also shown in Fig.l. It is an important feature of this approach 
that the calculation can be performed in the physical phase space as a large 
contribution of nonperturbative QCD arises from the extension of phase space 
from the quark level to the hadron level. It should be pointed out, however, 
that in theoretical calculations we take Mv = frir since we assume quark- 
hadron duality in our approach. 

Important properties of the distribution function are derived from field 
theory. Due to current conservation, it is exactly normalized to unity with 
a support < ^ < 1. It obeys positivity. When the distribution function 
becomes the delta function, — mft/M^), the free quark decay is repro- 
duced. Furthermore, the next two moments of the distribution function can 
be estimated in the HQET, as we shall discuss below. These two moments 
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determine the mean value /i and the variance of the distribution func- 
tion, which characterize the position of the maximum and the width of it, 
respectively: 

/x^Mi(0) = e + Mi(O, (18) 
^ M, ifi) = M2 (0 - (0 , (19) 

where M„(^) is the nth moment about a point C, of the distribution function 
defined by 

M„(o= (20) 

Jo 

By definition, Mo(0 = 1- 

The accuracy of the theory is remarkably improved by estimating the next 
two moments of the distribution function in the framework of the HQET. 
However, it cannot yet be completely determined in QCD. For practical cal- 
culations, therefore, we shall use an ansatz for the distribution function, 
which respects all known properties, with two parameters a and b as follows 

m = Nj^^^^^emi-o, (21) 

where N is the normalization constant. For a = rrih/MB and 6 = 0, eg . (121]) 
becomes a delta function, — rrib/MB), and the free quark decay is re- 
produced. Another form of the distribution function has been proposed in 

0- 



The perturbative QCD corrections to 0{as) has been calculated |jT5|, |TG 
It has the form 



2a 



^ rtp.rt „ H 



pert o 1 J 

on \mti J 
where is the free quark decay width: 



Mr. (22) 



r, = Fo$p , (23) 



with 

Glml IK 



$(x) = 1 - Sx^ + 8x^ - X® - 24:xHnx. (25) 
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An analytic expression for H{mc/mb) is given in |T6| . 

The semileptonic decay width can be calculated by substituting ([T0|) and 
( P^D into (^. The parameters involved and hence the sources of the theoret- 
ical error are: 

(1) the parameters in the distribution function (for the ansatz (|2l|) that will 
be used, they are a and b), 

(2) the beauty and charm quark pole masses nih and nic, 

(3) the strong coupling constant a^. 

There are several theoretical constraints on these parameters stemming 
from the HQET, which reduce the theoretical uncertainties considerably. We 
discuss them in turn. 

I. Performing a light-cone OPE and following the method of to expand the 
matrix elements of the local operators in the HQET, the next two moments 
and hence the mean value /i and the variance o"^ of the distribution function 
can be related |]1T| to two accessible parameters Ki, and Gb up to the order 

of {AqcD/mb)^, 

. = ^(l + B.). (26) 

where Eb = Kb + Gb- Both Kb and Gb are of order {AgcD/^b)'^- This leads to 
a model-independent conclusion: the distribution function is sharply peaked 
around ^ = /i ~ nib/ Mb and its width is of order Aqcd/Mb- 

For numerical analyses we need to know Kb and Gb quantitatively. The 
parameter Gb is related to the observables, 

m,Gfc = -^(MB. -Ms), (28) 

where the mass difference of the vector B* and the pseudoscalar B mesons 
is measured to be Mb* — Mb = 0.046 GeV. Kb can be reexpressed in terms 
of another often used parameter Ai instead. 

It is harder to determine Ai (or Kb). The accurate value of it is not known. 
Consequently, the mean value fi and the variance of the distribution func- 
tion are determined by the two parameters rrib and Ai: fi depends on rrib 
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strongly and Ai very weakly, while cr^ is sensitive essentially only to Ai. 
Hence, the parameters a and h in the ansatz (^) for the distribution func- 
tion are also determined by mf, and Ai. 

II. The quark mass difference is related to Ai in the HQET p] 

- = (Mb - Mz,)|i - + 0{l/ml)^, (30) 

where the spin-averaged meson masses 

Mb = ^{Mb + 3Mb*) = 5.31 GeV, (31) 
Md = ^{Md + 3Md') = 1.97 GeV (32) 

Finally, the remaining theoretical input parameters for our analysis are 
nib, Ai, and a^. 



3 Analysis 

We evaluate Tnonpert and Tpert in the approach described above using the 
three input parameters. For mb we use 

rub = 4.9 ± 0.2 GeV. (33) 



According to a QCD sum rule calculation |jT^, we take 

Ai = -(0.5±0.2) GeVl (34) 

As a result, the mean value and the variance of the distribution function are: 

/i = 0.93 ±0.04, (35) 
= 0.006 ±0.002. (36) 

A truncating of perturbative series causes the dependence of perturbative 
calculations on the renormalization scale fir- For inclusive semileptonic B 
decays perturbative QCD corrections are known only to the leading order. 
The result given in (p2D exhibits an implicit scale dependence of the strong 
coupling as. We vary the scale over the range of mb/2 < fir ^ f^b to estimate 
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the theoretical error due to the choice of the scale used in the argument of 
as. 

The dependence of the decay width F on the parameters is shown in Fig. 2. 
The variation of F with or Ai is stronger than /i^. The variation of mi, leads 
to an uncertainty of 8% in the decay width if other parameters are kept fixed. 
The same uncertainty in the decay width results from the variation of Ai. An 
uncertainty of 2% in the decay width is introduced when the renormalization 
scale fir is varied between mb/2 and rrib. In addition, the impact of the shape 
of the distribution function on the value of the decay width is studied. The 
value of the decay width is more sensitive to the variation of the mean value 
than the variation of the variance of the distribution function. Furthermore, 
we modify ( PT| ) with two more parameters a and (3 to be 

m = N^^^i^^^0K)0(l-O- (37) 

Using (37) we find that the value of the decay width is insensitive to the 
change of the shape of the distribution function if the mean value and the 
variance of it are kept fixed. This insensitivity diminishes the model depen- 
dence. This analysis yields 

7c = 49 ± 9 ps-i, (38) 

with a theoretical error of 18%. 

In Fig. 3, we compare the decay widths calculated in our approach, the 
free quark decay model, and the heavy quark expansion approach. The result 
in our approach shows that nonperturbative QCD contributions enhance the 
decay width by +(9 ± 6)% with respect to the free quark decay width, in 
contrast to the result of the heavy quark expansion approach where a reduc- 
tion of the free quark decay width by —(4.3 ± 0.5)% is found. The change 
of the sign indicates that the nonperturbative effects receive a large phase- 
space enhancement. We also observe that the decay width calculated in our 
approach goes to a free quark decay limit as — Ai decreases. This behavior is 
expected since the distribution function approaches a delta function, which 
reproduces the free quark decay, as — Ai and hence decrease. Thus this 
behavior provides a check of calculations. 

This theoretical analysis can be used to determine \Vcb\- Experimentally 
the inclusive semileptonic branching ratio Bsl has been measured at the 
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T(4S') and resonances, respectively. The lifetime tb has been measured 
by experiments at Z° and in pp collisions. The average of these measurements 
leads to |T^ 

= 67.3 ± 2.7 ns-^ (39) 

Putting it together with the theoretical value of 7c given in (|38|), we obtain 
from (I) 

\Vcb\ = 0.0371 ± 0.0007 ± 0.0034, (40) 
where the first error is experimental and the second theoretical. 



4 Conclusion and Discussion 

We have calculated the semileptonic decay width of the B meson using an 
approach based on the light-cone expansion and the HQET. Nonperturbative 
QCD effects are described by a single distribution function. Several impor- 
tant properties of the distribution function are known from QCD and the 
HQET of it. However, one still has to model the distribution function. For- 
tunately the result of the calculation of the decay width in this approach is 
nearly model- independent, since it is essentially only sensitive to the mean 
value and the variance of the distribution function, whose theoretical es- 
timates exist. Moreover, this approach is able to take into account both 
dynamic and kinematic components of nonperturbative QCD effects. We 
have shown that including the latter is indeed quantitatively crucial, which 
could increase the decay width significantly. We find an enhancement of the 
free quark decay width by +(9 ±6)% due to nonperturbative QCD contribu- 
tions, contrary to the claims from the heavy quark expansion approach. As 
a result, a value of \Vcb\ is extracted from the inclusive semileptonic B meson 
decay with a controlled theoretical error. 

The main theoretical uncertainty arises from the values of the b quark 
mass and the HQET parameter Ai. It seems possible to reduce theoretical 
uncertainties by a detailed fitQ to the measured charged-lepton energy spec- 
trum to determine the parameters and a calculation of the next-to-leading 
order perturbative QCD correction.|^ Future measurements of the distribu- 

"''Such a fit has been done in the heavy quark expansion approach ||l9| . 
^Partial calculations of higher-order corrections exist po|. 
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tion function and more theoretical efforts on calculations of hadronic matrix 
elements should enable to further reduce the uncertainties. 

Careful inclusion of the kinematic effect of nonperturbative strong inter- 
actions is also necessary for reliable predictions for the nonleptonic decay 
widths of hadrons containing a b quark. The nonleptonic decay widths of b 
hadrons may be calculated in a similar way. The decay widths are expressed 
in terms of the b hadron masses rather than the b quark mass provided the 
phase-space effects are included, whereas according to the heavy quark expan- 
sion, the relevant mass in the decay widths should be the universal 6-quark 
mass and no corrections of order l/rrih should be present We would an- 
ticipate an enhancement of the nonleptonic decay width by nonperturbative 
QCD if both dynamic and kinematic effects of it are properly taken into ac- 
count. Since the phase-space effects cancel out to a large extent in the ratio 
of the decay widths, they are unlikely to significantly change the semileptonic 
branching ratio of the B meson or the average number of charmed hadrons 
produced per B decay. On the other hand, the prediction on the ratio of the 
Afc and B lifetimes from the heavy quark expansion approach seems to be in 



conflict with the data ||2^. In this case phase-space effects are enlarged due 
to the significant difference between the A;, and B masses. The replacement 
of the 6-quark mass with the non-universal 6-hadron masses results in a per- 
fect agreement between the theory and the experimental data, giving 
evidence for the phase-space effects. 
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Figure 1: Phase space for the 6 — > c inclusive semileptonic decay. The interior 

of the sohd curve is the hadron level phase space (the changeablcncss of the 
mass of the final hadronic state is not shown explicitly). The interior of the 
dashed curve is the quark level phase space. 
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~Ai (GeV^) 

Figure 2: Dependence of the semileptonic decay width F on the theoretical 
input parameters rrib, Xi, and Hr- The sohd (dashed) curves are for the 
renornialization scale /ir = {l^r = irih/'^)- The curves with solid dots, 
boxes, triangles correspond to rrib — 4.7, 4.9, 5.1 GeV, respectively. 



16 



0.2 0.3 0.4 0.5 0.6 0.7 

Figure 3: Semileptonic decay width F as a function of Ai calculated in our 
approach (solid curve), the free quark decay model (dotted curve), and the 
heavy quark expansion approach (dashed curve). We take = 4.9 GeV 
and as — 0. 
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